AD-784  894 


SCHEDULING  INDEPENDENT  TASKS  ON 
NON -IDENTICAL  PARALLEL  MACHINES 
TO  MINIMIZE  MEAN  FLOW- TIME 

Douglas  C  lark 

Ca  rnegie  -  Me  1  Ion  University 


Prepared  for: 

Defense  Advanced  Research  Projects  Agency 
Air  Force  Office  of  Scientific  Research 


June  1974 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Sprinpfield  Va.  22151 


y 


UNCLASSIFIED  _ 

JWITV  CLaSS  FiCATiqh  This  PAGE  f»hin  lie,. 

REPORT  DOCUMENTATION  PAGE 


1 .  hi  PORT  HU  U  fl  I .  f  I 


AfOSR  -  !►;  -74  -  1  3  5  % 


PAGE  KM  AD  INSTRUCTIONS 

- _ nF.ro;;fC  COMi'U  "INO  FORM 

P  OOVT  *CCSilCh  MO  i  RECU  icn  vscaTZlog  number - 


*•  Tl  TLC  (t»nd  Subtitle) 

SCHEDULING  INDEPENDENT  TASKS  ON  Nl.vl-1  DENT! CAL 
TARALLEL  MACHINES  TO  MINIMIZE  ML.A'I  FLOL'-TIHt 


|7  AUTHOR'* 


Douglas  Clark 


1  TYPE  or  REPORT  «  PERIOD  COVfNEO  j 

Interim  J 

t  PERFORMING  ORG  REPORT  NJMHER 
6.  CONTRACT  OR  GRANT  NUMBER!*) 

F44620-73-C-0074 


9  performing  organ  7 at  ioi,  NAyr  and  address  ~<^u7trvUT^7T,. NAuTT-T7sT 

Co rneg ic-Mc  ]  Ion  Vni  wr^ity  ar^a  6  work  umt  numbfrs  ’ 

Pcpar traer.L  of  Computer  Science  61101D 

Pittsburgh,  PA  15213  A0Z466 

"  C,0NT”0LV'!C  ^  MCf  NAMF  *N0  ADDf,C“  ‘  '  "  Tt  RLPORT'TTTi - 

Defense  Advanced  Research  Projects  Agency  Dine  1974 

1400  Wilson  Blvd  - * - 

Arlington,  VA  22209  ^  9 

'  Air  ForcoC  Ot  f  icV'oi"  Scienti  fie" "KarclT"'0';'' «cu,.,t7cua5s.  m 

1400  Wilson  Blvd  /  UNCLASSIFIED 

Arlington,  V'A  22209  _ _ 


*1. 


UNCLASSIFIED 


—— —  _  _  ’ 

■So.  I>[. CLASSIFICATION  'DOWNGRADING  J 

SCHEDJLE  I 


n,  oisTRifu’-io',;  sv  M 1  silnt  <o i  mi.  ^ , 


Approved  for  public  release;  distribution  unlimited. 


17.  ni?TR,0uTlON  STATEMENT  7c/  Ibe  abstract  entered  In  lilork  30.  II  dlller.nl  fr,.ri  Heport) 


ic.  supplementary  notes 


IS  KEY  WORDS  fCo-  'Inna  cn  r.-.oi..  .10,  II  nacaaamy  mcf  identity  Ty  blo.k  nunJZ) - - - - - 

Reproduced  ty 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  S  Dopartmenl  of  Commerce 

To — IbLto.  T  — : . . . . _  Sp'mKfield  VA  22151 

A  r  ill  1  r7-  h  -i  ^  ed*  II  n*ce»*ery  and  Identity  by  block  number ) 

lr  i  f  tasks  having  known  processing  time  requirements  on  a  set  of 

n;  £;  r1  ParallC  rhineS  iS  tC  be  Scheduled  that  the  mean  flow-time 
of  tnc  tasks  is  as  small  as  possible.  In  this  paper  it  is  shown  that  a  tri- 

*  extension  of  a  simple  algorithm  for  a  restricted  case  perforas  well  and 
often  epemany,  in  the  general  case.  A  principal  result  is  that  for  every 
problem,  some  renumber mg  of  the  tasks  will  cause  this  algorithm  to  produce 

rithmVr?1  SCb0dule:  L'[’per  bounds  on  the  worst-case  performance  of  the  alp0- 
_nthm  „re  gi.cn,  and  average  performance  is  explored  u'siag  Monte  Carlo  tcchniqi  es 

UD  I  jan  71  1473  LDi  TiO'  or  I  NOV  es  IS  OBSCLCTE  ' 

UNCLASSI1  TED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  /Wh  enhata  tnlerrj’ 


II 


SCHEDULING  INDEPENDENT  TASKS  ON  NON- IDENTICAL 
PARALLEL  MACHINES  TO  MINIMIZE  MEAN  FLOW-TIME 


Douglas  Clark 


Department  of  Computer  Science 
Carnegie-Mellon  University 
Pittsburgh,  Pa.  15213 

Jur?e  1974 


This  work  was  supported  in  part  by  a  grant  from  the  Xerox  Corporation 
Palo  Alto  Research  Center  and  in  part  by  the  Ad  'anced  Research  Projects 
Agency  of  the  Office  of  the  Secretary  of  Defense,  Contract  Number 
F44620-73-C-0074,  which  is  monitored  by  the  Air  Force  Office  of  Scien¬ 
tific  Research.  This  document  has  been  approved  for  public  release  and 
sale;  its  distribution  is  unlimited. 


ABSTRACT 


A  collection  of  tasks  having  known  processing  time  requirements  on  a  set  of  non-ldentlcal  parallel  mrchlnen 
Is  to  be  scheduled  so  that  the  mean  tow-time  of  the  tasks  Is  as  small  as  possible.  In  this  paper  It  Is  shown 
that  a  trivial  extension  of  a  simple  algorithm  lor  a  restricted  case  performs  well,  and  often  optimally,  in  the 
general  case.  A  principal  result  Is  that  for  every  problem,  some  renumbering  of  the  tasks  will  cause  this  al¬ 
gorithm  to  produ.e  an  optimal  schedule.  Upper  bounds  on  the  worst-case  performance  of  the  algorithm  are  given, 
and  average  performance  Is  explored  using  Monte  Carlo  techniques. 

1.  INTRODUCTION 

Tlie  problem  addressed  In  this  paper  Is  the  sequencing  of  n  Independent  tasks  on  m  parallel  and  non-ldentlcal 
machines  so  that  the  average  flow-time  of  the  tasks  is  as  small  as  possible  [5],  We  will  assume  that  all  n  tasks 
or  lobs  are  simultaneously  available  at  time  rero;  that  there  are  no  feasibility  or  precedence  constraints  among 
the  tasks;  that  tasks  may  not  be  preempted;  that  a  machine  can  process  only  one  job  at  a  time;  and  that  the  pro- 
cessti  q  time  required  by  a  Job  1  on  a  machine  J  Is  given  by  a  positive  number  p^.  The  Inability  of  some  machine 
to  process  some  task  may  be  represented  by  making  the  corresponding  p^  prohibitively  large. 

Fig.  1(a)  Is  a  processing  time  array  p  for  an  eight  Job,  three  machine  problem.  A  schedule  for  thla  problem 
Is  shown  In  Fig.  1(b)  In  the  form  of  a  Cantt  chart  [A],  which  Illustrates  the  p'.-allel  activity  of  the  three  ma¬ 
chines  along  a  horiroi. tal  time  axis.  The  rectangular  blocks  in  the  chart  have  len.ths  equal  to  the  processing 
ilmea  of  the  Job.  with  whose  numbers  Ihey  are  labeled.  The  flow-time  or  t lme-tn-sy  Uem  of  a  Job  in  a  particular 
schedule  Is  simply  the  time  at  which  that  Job  completes  its  execution,  where  the  schedule  begins  at  time  0.  Thus, 
In  Fig.  1(b)  the  flow-time  of  Job  3  Is  2,  of  Job  A  is  10,  of  Job  6  Is  12,  and  so  on. 

Borrowing  some  notation  from  [5],  We  w  11  denote  by  Jf 1 ]  the  number  of  the  Job  scheduled  1th  on  machine  J. 

In  Hg.  1(b),  1 T 1  1  1»  5,  1  r2  ]  is  A,  1[3]  is  2,  and  so  on.  Let  fj  be  the  flow-time  of  Job  1  In  »  particular 
schedule,  let  F  be  the  sum  of  the  flow-times  of  all  n  Jobs,  and  let  be  the  number  of  Jobs  scheduled  on  machine 
J.  It  Is  clear  that 

fJPJ  ’  PJ[  1  ] .  J 

fjp] "  pjn  j,j  4  pj[2i,j 


*nd  In  general , 


'  Pjm.J  4  P)[2],j  4  4  PJ[nJl,j’ 


fjm  ’  kt,  pj[ki.j- 


In  Up.  l(bj,  lor  example,  we  have 


1 


f1[l]  '  P1f1),l 


fH2]  '  pi;i],i  +  P1[2],1 


f'm  '  pl|l),1  +  P1  [2  ] ,  1  +  P I  [  3  ] ,  1  -  >  +  9  +  3. 


Hie  total  flow-time  of  a  schedule,  F,  may  be  expressed  as  a  sum  of  contributions  from  each  machine.  Collect¬ 
ing  terms  in  equations  (1),  we  may  express  the  contribution  of  machine  J  as 

n .  n 

y  yJ  f 

i-i  fJC‘l  '  tt,  k_, 


nJPJ[l).J  +  <nj"')pjr2],  '  +  ”•  +  2pJ[nj-l],j  +  pj[nJ],r 


Summing  over  all  machines,  we  get 


y  yJ 

T  m  L.  L  f  »  i_,  u,  p 
J  ’  J-l  i-1  k-1 


11 

m  j  i 

I  I  I 


J-!  i-1 


n1P1[l],1  +  (V,)pl[2],l  +  ”•  +  2pl[n,-1],l  +  P1  [n,  ] ,  1 

|  n2P2[l],2  +  (n2-,)p2[2].2  +  "•  +  P2[n2],2  + 

+  nmPm[1],m  +  ( V 1  )Pra[2  ],m  +  +  Pm[n  j  „• 

m 

For  example,  in  Fig.  1(b)  we  have 


(2) 


F  -  3-1  +  2-9  +  1.3 
+  2-2  +  1.1 
+  3.3  +  2-4  +  1.5 
-  51. 

Our  goal  is  to  minimize  mean  flow-time,  “F,  but  minimizing  F  itself  is  equivalent,  and  we  will  follow  that 
approach  in  the  rest  of  this  paper.  Fig.  1(c)  gives  an  optimal  schedule  (optimal  schedules  are  not  necessarily 
unique)  for  the  problem  of  Fig.  1(a),  with  F  -  34. 

A  non-enumerative  algorithm  or  minimizing  mean  flow-time  in  the  general  case  was  discovered  by  Bruno, 
toffman,  and  Sethi  2,3].  Their  algorithm  is  based  on  a  reduction  of  the  problem  to  a  minimum-cost  network  flow 
problem,  and  the  time  required  by  the  algorithm  is  0fn3)  when  n  2  m  (the  case  of  interest)  and  0(n2m)  when  n  <  m. 
This  paper  proposes  and  analyzes  an  algorithm  which  finds  schedules  that  are  good,  and  frequently  optimal,  with 
respect  to  mean  flow-time,  and  does  so  at  very  small  computational  co9t. 

Section  2  of  this  paper  reviews  an  easy  algorithm  for  an  important  restriction  of  the  general  problem.  In 
section  3  this  algorithm  is  extended  to  cover  the  general  case,  and  it  is  shown  that  while  an  optimal  schedule  is 
not  always  produced,  the  performance  of  the  algorithm  strongly  depends  on  the  ordering  of  the  rows  of  the  proces¬ 
sing  time  array  P.  Section  4  examines  analytically  the  worst-case  performance  of  the  algorithm  under  various 
ordering  rules,  and  the  algorithm’s  average  performance  under  these  rules  is  explored  empirically  in  Section  5. 
Section  6  contains  the  conclusions  of  the  paper. 
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2.  MACHINE  FACTOR  CASE 

An  Important  and  realistic  restriction  of  the  general  problem  erases  when  each  p^  is  the  product  of  a  time 
aaaoclated  with  job  i  and  an  efficiency  fa(  tor  easociated  with  machine  J,  that  is,  p{  (  -  p^.  In  this  restrict¬ 
ed  case,  mean  flow-time  is  minimised  by  a  simple  procedure  which  follows  immedietely  fro™  the  analysis  of  [5] 
First,  rewrite  equation  (2)  as  follows: 


m  J  i  m  J  1 

F  -  I  I  Ip  .HI., 

j-1  i-1  k-1  j.]  la1  k»,  i  J[kl 


'  niwipu i ] +  +  -  +  2-,p, +  vun,] 

+  n2W2p2[1]  +  <nr,)w2P2[2]  +  +  W2P2l  u_]  + 


+tln>w.pmrn  +  (n-.-,>*Lp  w  p  . 

tnimLlJ  m  mm2]  mKm[n  ] 

m 

Second,  since  F  is  a  sum  of  n  terms,  each  of  which  is  e  product  of  one  of  the  n  p^s  with  one  of  the  n  coeffici¬ 
ents  n]wl ,  (n]-1)w1 . w],n2W2 . w2 . ,nmWm,“-,Wm’  pick  e*  coefflci««ts  the  n  smallest  of  the  nm  possibil¬ 
ities  nWj ,  (n-l)Wj . w)  >nw2 . W2  * . "“m . Wm'  Thi*  wil1  determine  the  values  of  the  nJt  Third,  mini¬ 

mise  F  by  metching  the  lergeat  ?1  with  the  smallest  coefficient,  the  second- lerge.t  Pi  with  the  second-smallest 
coefficient,  end  so  on.  If  e  particular  pt  is  matched  with  (n^-klw^,  then  job  i  is  scheduled  (k+l)th  on  machine 

J.  Picking  coefficien*  (n^-kjw^  implies  that  we  have  already  picked  (n^-k-Dw^ , (nJ-k-2)wJ . 2Wj,  end  Wj,  ao 

the  resulting  schedule  is  well-formed.  (Tiat  is,  it  is  not  possible,  for  example,  to  schedule  some  Job  fourth 

on  some  mechine  end  not  schedule  some  oth.,r  jobs  first,  second,  f  nd  third.) 

Eech  choice  of  a  coefficient  can  bj  res,*-lcted  to  be  emong  only  m  possibilities  out  of  the  totel  of  nm.  The 
first  chosen  will  surely  be  on.  of  w, . «m  uhe  smell.st,  in  fact).  Suppose  it  is  w^  Then  the  second  coeffici¬ 
ent  will  be  the  amellest  of  w, . 2V"’V  At  each  stage  the  integer  multiplier  of  the  chosen  coefficient  is 

increesed  by  1,  end  the  next  choice  made.  The  schedule  is  thus  being  determined  "back"  to  "front". 

An  exemple  of  this  procedure  is  given  in  Fig.  2(e)  for  a  five  job,  three  mechine  problem.  The  Jobs  happen 

to  be  numbered  so  thet  P,  *  P2  *  ...  *  P5 ;  consequently  no  sorting  of  the  pt  is  necessary.  In  eech  row  of  the 

teble  the  smallest  of  the  m  pot.ntiel  coefficients  is  chosen  (circled),  and  in  the  next  row  the  circled  coeffici¬ 
ent  is  increesed  by  the  corresponding  Wj.  In  row  2  there  ere  two  smellest  coefficients;  the  choice  bnween  them 
is  arbitrer  .  since  different  schedules  with  the  same  (optimal)  F  will  result  from  different  choices  in  the  case 
of  a  tie.  Fig.  2(b)  shows  the  resulting  schedule  and  the  celculatlon  of  F. 

Fig.  2(c)  shows  the  application  of  en  equivelent  algorithm  to  the  problem  of  Fig.  2(a).  Here  the  processing- 
time  errey  P  is  explicitly  shown,  end  the  coefficient,  in  each  row  now  represent  the  possible  sequence  positions 
on  the  three  mechines,  counting  from  the  end  of  the  sch.  ,’ule.  for  the  corresponding  Job  It  is  important  to  n.  te 

thet  for  each  column  J,  P]J  i  p  r  ...  *  p  Let  the  coefficients  in  e  particuler  row  i  be  h.,h, h  (ir 

this  cese,  m  -  3).  Then  the  minimum  h^p^  is  chosen,  ptJ  is  circled,  and  hj  is  increesed  by  1  in  the  next  row. 


machine  factors: 


2.5,  w3  .  1 


V,  - 


2, 


w„  ■ 
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This  Is  exactly  equivalent  to  the  algorithm  cf  Fig.  2(a),  and  the  same  schedule  results.  The  difference  lies  in 
the  lack  of  explicit  use  of  the  w^  in  the  second  version  of  the  algorithm,  and  this  difference  will  be  exploited 
in  Section  3. 

Sorting  the  p^  requires  computational  exertion  of  0(,n  log  n) ,  and  finding  the  smallest  coefficients  requires 
0<;nm),  so  the  time  complexity  of  the  machine  factor  algorithm  is  0(max(n  log  ..,nm)). 


3.  EXTENSION  TO  THE  GENERAL  CASE 

The  second  version  of  the  procedure  of  the  preceding  section  (Fig.  2(c)),  since  it  does  not  use  the  machine 
factors  w^  may  be  applied  in  the  general  case,  though  with  no  guarantee  that  optimal  schedules  will  be  found. 

Fig.  3(a)  shows  the  result  of  applying  this  procedure  to  the  problem  of  Fig.  1(a).  Since  this  problem  is  outside 
the  machine  factor  case  we  cannot  sort  the  jobs  according  to  a  single  processing  time,  so  in  Fig.  3(a)  the  job 
numbering  of  Fig.  1(a)  is  retained,  and  the  schedule  of  Fig.  3(b)  results.  The  value  of  F  is  not  optimal,  but 
the  schedule  was  very  quickly  computed,  so  this  algorithm  will  hereafter  be  called  the  Quick  And  Dirty  (QAD)  al¬ 
gorithm.  The  coefficients  used  at  a  particular  stage  of  the  schedule  generation  will  be  called  the  QAD  coeffici¬ 
ents  and  labeled  h. ,h, . h^.  In  each  row  i  of  P  the  QAD  algorithm  chooses  the  minimum  of  h jPjq • h2P i 2  *  * ' * ,hmPim 

and  Increments  the  chosen  h^  in  the  following  row.  As  in  the  machine  factor  case,  the  choice  of  hjPtj  means  that 
job  1  is  scheduled  (n^-h^+llt^  (or  equivalently,  h^th  from  the  end)  on  machine  j,  and  tha-  the  term  h^p^  will 
occur  in  the  calculation  of  F.  (This  latter  fact  is  very  Important  for  the  rest  of  this  paper.) 

An  obvious  infirmity  of  the  schedule  in  Fig.  3(b)  is  that  the  Jobs  are  not  scheduled  in  "shortest  processing 
time  first"  (SPT)  order  on  each  machine.  SPT  order  is  optimal  for  the  mini'* 'ration  of  F  on  a  single  machine  [5], 
so  an  optimal  schedule  for  the  more  general  problem  clearly  must  have  Jobs  in  SPT  order  on  each  individual  machine. 
This  fact  suggests  that  the  performance  of  QAD  may  be  improved  by  adding  to  it  a  procedure  which  sorts  the  jobs 
in  SPT  order  on  each  machine  after  the  schedule  has  been  produced.  Call  the  combined  procedure  QAD*.  If  QAD*  is 
applied  to  Fig.  3(a),  the  result  is  a  schedule  with  an  improved  F  of  35,  which,  while  still  not  optimal,  is  con¬ 
siderably  better  than  QAD's  F  ■  48. 

Though  an  optimal  schedule  was  not  found  in  Fig.  3(a),  if  the  rows  of  P  are  permuted  as  in  Fig,  3(c),  the 
optimal  schedule  of  Fig.  1(c)  will  be  generated  by  QAD  (and  therefore  also  by  QAD*).  This  might  appear  at  first 
glance  to  be  no  more  than  a  fortunate  coincidence;  the  following  theorem,  however,  states  otherwise. 


Theorem  3-1 .  Given  an  arbitrary  processing  time  array  P,  there  exists  a  permutation  of  the  rows  (a  renumbering 
of  the  jobs)  such  that  the  QAD  algorithm,  operating  on  the  permuted  array,  will  yield  a  schedule  with  optimal 
(minimal)  F. 


Proof.  Let  P  be  arbitrary.  The  plan  of  the  proof  will  be  to  show  that  if  the  first  k  jobs  (Jobs  l,2,...,k)  of  P 
are  scheduled  optimally  by  QAD,  then  one  of  the  remaining  n-k  jobs  can  be  renumbered  k+1  so  that  it,  too,  will  be 
scheduled  optimally.  If  this  can  be  demonstrated  for  k  ■  0, 1 ,2 , . . . ,n-l ,  then  the  theorem  will  be  proved.  In 
this  proof  the  phrase  "scheduled  optimally"  will  mean  scheduled  (by  QAD)  i_n  agreement  with  some  optimal  schedule . 
both  in  the  assignment  of  jobs  to  machines,  and  in  the  order  in  which  they  are  assigned.  Recall  that  QAD  assigns 
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Jobi  from  last  to  flrat  on  each  machine. 

Let  S  be  an  optimal  achedule  for  P.  Aaaume  that  the  flrat  It  Jobe  of  P  are  scheduled  by  QAD  In  agreement 
with  S,  where  k  may  be  any  Integer  from  0  to  n-1.  (There  la  nothing  to  prove  If  k  ■  n.)  We  will  try  to  find  a 
(k+l)th  Job  such  that  either: 


1,  Jobs  1  through  k+1  will  be  acheduled  by  QAD  In  agreement  with  S;  or 

2.  there  la  some  other  optimal  schedule  S'  such  that  Jobs  1  through  k+1  will  be  scheduled  In  agreement 
with  It, 


If  at  least  one  of  these  two  alternatives  la  always  true,  the  theorem  Is  proved. 

Fig.  4  Illustrates  the  situation  with  an  example.  The  shaded  jobs  are  Jobs  1  through  k,  optimally  scheduled 
by  QAD.  If  the  QAD  schedule  Is  to  agree  with  S,  the  only  Jobs  which  are  "candidates"  for  Job  k+1  are  those 
marked  by  an  asterisk  in  Fig.  4.  There  Is  at  most  one  such  candidate  Job  on  each  machine.  Let  g^  (J  "  1,2 . m) 


be  the  number  of  the  candidate  Job  on  machine  J  if  there  Is  one;  otherwise  let  g^  •  0. 
the  QAD  coefficient  for  machine  J  in  row  k+1  of  P. 


Let  hj  (J  ■  1,2 . m)  be 


Figure  4 


Cons 


truct  a  graph  with  m  vertices  v^.v^ . v^.  For  each  vertex  draw  a  directed  arc  (vt»vj)  from 


to 


v^  if  and  only  if 


1.  gt  t  0;  and 

2.  QAJ  would  schedule  job  gt  on  machine  J  If  It  became  Job  k+1  of  P. 


The  existence  of  an  arc  (v  ,v  )  implies  that  h  p  s  h  p  ,  with  equality  only  if  the  tie-breaking  rule  used 

*  J  J  o^iJ  1 

by  the  QAD  algorithm  would  choose  machine  J  over  machine  1  for  Job  g  . 

If  there  is  any  arc  of  the  form  (v1>vt)i  then  Job  gt  can  become  Job  k+1,  QAD  will  schedule  it  according  to 
S,  and  we  are  done.  Assume,  therefore,  that  there  are  no  such  arcs  in  the  graph. 

Suppose  some  vertex  v^  has  one  or  more  arriving  arcs  but  no  departing  arc,  and  let  Vj  be  a  predecessor  of 
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vj  (^8-  5(b)).  Construct  a  new  schedule  S'  which  Is  Identical  to  S  except  that  Job  g^  Is  scheduled  first  on 
machine  J  (Fig.  5(a)).  If  we  constructed  the  graph  corresponding  to  S’,  we  would  find  Fig.  5(c)  In  place  of 
Fig.  5(b).  If  S'  Is  optimal,  then  Job  gf  could  become  the  new  Job  k+l  In  P  and  be  optimally  scheduled. 


Figure  5 


Denote  by  F 


S 


and  F 


S' 


the  total  flow-times  of  the  two  schedules. 


Then  we  have 


„S '  8 

F  -  F  +  h  p  .  -  h  p  -  c. 

J  Sj.J  1  gj.l 


(3) 


Hie  non-negative  term  c  occurs  in  (3)  because  the  removal  of  Job  gt  from  machine  i  cause t  the  coefficients  of  all 
preceding  Jobs  (If  there  are  any)  to  be  reduced  by  1  (see  Fig.  5(a)).  The  existence  of  arc  (v^v^)  Implies  that 
hJPg,.J  S  hiY,l*  ,nd  the  optimality  of  S  Implies  that  FE  2  FS.  We  conclude  from  (3),  therefore,  that  c  -  0, 

s'  1  s  1 

F  *  F  ,  and  S'  is  optimal.  Thus  Job  can  become  Job  k+1  of  P. 


Now  suppose  that  the  graph  contains  no  vertices  that  have  arriving  arcs  but  no  departing  ere.  Since  there 
Is  at  least  one  arc  lr.  the  graph  (becaune  k  <  n)  ,  since  the  number  of  vertices  Is  flnlti,  and  since  there  are  no 
arcs  (v^),  we  conclude  that  there  must  exist  a  cycle  In  the  graph.  Suppose  the  cycle  has  three  arcs,  (v^Vj) 
(vJ’V’  *nd  (vk’vl)*  **  *hown  ln  Fi8-  6<b>-  (Hie  following  argument  may  easily  be  generalized  to  cycles  of  any 
size.)  Construct  a  new  schedule  S'  Identical  to  S  except  that  Job  gt  moves  to  Job  g^ ’ s  position  on  machine  J, 
Job  gj  goes  to  machine  k,  and  Job  to  machine  1  (Fig.  6(a)).  n,e  graph  corresponding  to  S'  would  Include  Fig. 
6(c)  in  place  of  Fig.  6(b). 

S  S 1 

With  F  and  F  as  before  we  may  write 


S 1  S 

F  -  F  +  (h.p 


j V  '  *  lvV '  "‘V1* 4  '"‘V1 '  vtk.k>- 


(4) 


The  existence  of  the  arcs  ln  the  cycle  means  that  each  parenthesised  term  ln  (4)  is  non-positive,  and  since 
S'  S 

F  2  F  ,  we  conclude,  as  before,  that  S'  Is  optimal.  Any  one  of  Jobs  gj ,  and  gk  can  become  Job  k+l  in  P 
be  optimally  scheduled. 


and 


We  have  therefore  shown  what  we  set  out  to  show,  namely,  that  there  must  always  exist,  for  k  -  0,1 . n-1. 
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(b) 


(c) 


(•) 

Figure  6 

“  (k+1) th  job  that  will  be  optimally  acheduled  bv  OAD  Th.  .......  , 

b/  QAD  The  existence  of  a  row-permutation  of  P  that  yield,  an 

optimal  achedule  follow.  Immediately.  ■ 

One  consequence  of  Theorem  J-l  i.  that  an  optimal  .chedule  may  be  found  by  applying  QAD  to  each  of  the  n' 

'“T'"10"  ■' ' *■  *“*  -  - — -  - .........  .. 

’■  -“h  -  «”  or  „„ 

"“,4  "*  ™»  -<  P  W.  .o  .h.  execution  .1  w  „  qA0. 

““  ‘h*  .  «  «  -ruu.  ..  1..., 

'•  — *  1.  •  .eh~  tor  tlnOlu,  .1^1.  no...,  Pl  i.  ,h.  ^  t<pt„.nt> 

"  *  )0b’  "*  Vh‘"'  "•  U“''  “  *  k«P  '»  .f  P.  Scv.ral  .u.h  .or,,., 

rule,  are  proposed  and  analysed  1„  the  next  two  .ectlon.  4f  thl.  paper. 

4-  WORST  CASE  PERFORMANCE 

section  we  will  examine  the  worst  case  performance  of  the  QAD  algorithm  under  various  row-sorting 
mi...  Sloe.  QAP*  ..K...1.  m  ^  ^  b>  ^  ^  >ni 


“•  ‘““'P  'hith  *"  «■»  <"  W  ■!.  not  nec....rlly  „ 

"  “*  '“'1  W  '»  •“*  '  r».  .,.  „  „„ 

*“■  “  P“"  “  ““  “*“>  “”-“™  -  '•  -  — »  •«  ....  .  .in  ,. 


RULE 

QAD  . 

F 

OPT 

-  -Ill  «r«  ....  uPP.,  OounO.  „  ,xPr...,„n  f„  various  s>  ..  ^ 

Crlo  technique,  to  estimate  the  average  performance  of  th.  algorithms. 

Ut  P1  .„d  P2  be  two  row-permut.tlon.  of  an  arbitrary  P.  Denote  by  F1  and  F2  the  total  flow-time, 
of  their  QAD  scedules.  Then 

f’ 

~r  <  n 
F 


(5) 


And  this  Is  a  best  bound. 


12  12 
Proof.  Let  and  denote  the  contributions  of  Job  1  to  F  and  F  ,  respectively,  so  thst 


C1  *  hJPiJ 
c\  “  hkptk 


(6) 


1  2  1  ? 
for  some  j  and  k,  where  end  ^  are  the  chosen  QAD  coefficients  for  Job  1  in  P  and  P  .  We  will  show  that  for 

all  1, 

.1 


-j  s  n. 
C1 


(7) 


Suppose  (7)  is  false  for  some  i.  Then,  using  (6),  we  may  write 


h]pu >  ^ik- 


(8) 


By  the  operation  of  the  QAD  algorithm,  we  know  that 


and  In  particular. 


h|p  -  min  (h 1 p  ) 
i  1J  ISgft,  8  ig 


h>iJ  ShkPlk 


(9) 


where  k  la  as  In  (6).  Combining  Inequalities  (I  )  and  '#),  we  get 

^  hlPlJ  >  "'W 

hkPlk  >  "Vlk 
>  • 


(10) 


1  2 

Inequality  (10)  la  plainly  Impossible,  since  and  hj^  are  Integers  between  1  and  n.  This  contradiction  proves 


(7). 


Now  write 


Fi  c,1+c’+...+c1 

F  m  1  2  n 

F  2  2  2  * 

2  C.+C,+...+C 

1  2  n 


From  Inequality  (7)  and  the  fact  that  all  cj  and  are  positive.  It  follows  thst 


c]+c’+...+c1 

2 _ n 

"5  2  ? 

C.+C.+...+C 
1  2  n 


£  n. 


(11) 


For  equality  to  hold  in  (11),  It  must  be  true  that  for  all  i,  cj  -  nc|.  But  cler.ly  this  cannot  be,  so  Inequal¬ 
ity  (11)  becomes  strict  and  (5)  Is  proved. 

To  show  that  n  Is  a  best  bound,  let  P  take  the  following  form: 
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n-1 

n 


12  3 


*0)0) 
«  0)  0) 


0) 

0) 


02) 


I  It  ID  . 

X  (u  o>  • 


Ui 

0) 


where  «  <  <  X  <  <  u>.  Let  P  be  ae  given  and  construct  P2  by  interchenging  rows  1  and  n.  The  QAD  elgorlthm  will 
schedule  .11  Job.  on  machine  .  for  both  P>  and  P2;  only  the  Job  order  will  be  different.  He  will  have 


and 


L.  m  '*  +  +  3t  +  ...  +  (n-1 ) c  +  nX 
p2  X  +  2*  +  3i  +  ...  +  (n-1)*  +  n* 


.  4  F1  tC. 

lim  n. 

«-Kl  F  * 


Thu.  we  may,  for  m  *  1  and  n  2  2 ,  con.truct  en  errey  P  with  two  row-permut.tlon.  p’  and  P2  such  thet  the  ratio 


f‘/f  1.  arbitrarily  clo.e  to  n.  Therefore  n  1*  e  best  bound. 


CoroLlery  4-1.  Let  F^  be  the  total  flow-time  of  the  QAD  schedule  for  an  erbitrery  P  under  en  arbitrery  row- 
permute tion.  Then 


^<n 


*nd  this  is  a  beat  bound. 


OPT 


03) 


fraof.  Thet  n  1.  e„  upper  bound  follow,  directly  from  Theorem  3-1  end  Theorem  4-1.  That  it  1.  e  beat  bound  may 

be  won  by  examining  (12)  in  the  proof  of  Theorem  4-1.  In  that  example  F2  -  F  iJ 

OPT*  m 

'"*  b~«"  l"  *»«-  M  Coroll. ry  A-l  d.p.,d.  on  a,,  yooull.r 
tor.  of  th.  proce..tng>tlme  .rr.y  (11).  If  QAn*  „r.  u..d  ln.„.d  of  gu.  th.  .„i„.i  ,ch.d„l.  (,1)_SPT  „ 
ooehlno  .—Id  b,  fourd.  th.  boood  „  l„  (5)  .„d  (13)  tl  b„nd  „„ 

Tb.  ..o.lbl.  »„  of  both  elgorlthm.  d.p.„d.  o,  ......  ...  roo.  of  P  prior  to  tb.I, 

tioo.  Ib,  approach  ».  t.U  .1,1  b.  to  .or,  ,h.  ro..  .o  tho,  ,h.  o,  .cm.  .logu-tad  f„„.,,o„  of  ...h 

ro.  ar.  I.  non-lncr...,.,  ord.r,  P.rh„.  ,h.  mo.t  obvlou.  th.lt.  ,  tat,,.„  tt.  tl.  . 

proc...tng  .f  ,  job.  TMc  oth.r  ,o..lblll,l..  ...  ,b.  ..«lm„.  ,„d  .1„1.„.  ,p,  .  ,„y 

of  ,h...  thro.  ml..  .Ill  .....  W  perform  .,ch„.  . . .  1,  ..... 

require  thet  this  be  true  of  whetever  sorting  rule  we  edopt. 

coot,  by  Aim.  MAX,  .„d  MIN,  th.  ru...  .ortlog  lo  ord.r  of  «oh-l,cr.,.l„,  ro.  ...r.,., 

r..p.t„o.,y.  „.  ...  ...fly  .ho.,  unfortunately,  ,b.,  A»E  .,d  MAX  do  no.  ..pro,.  ,b.  ucr.t  ....  „tfor- 

mence  of  QAD.  For  consider  the  following  erray  P; 


13- 


fl 


i 

2  . 

a  •  m 

1 

c 

iv+X  . 

.  .  ui+X 

2 

c 

ur*X  . 

.  .  ui+X 

• 

• 

• 

• 

n-1 

< 

ur*X  . 

.  .  u)+X 

n 

X 

ID 

.  .  it 

where  t  <  <  X  <  <  id.  The  rowe  are  already  arranged  according  to  both  AVE  and  MA> .  The  optimal  achedule  for  (14) 
is  SPT  on  machine  1,  but  QAD  schedules  Job  n  first  instead  of  last  on  machine  1.  Therefore,  using  the  argument 
from  the  proof  of  Theorem  4-1,  the  bounds 


_avf 

5^<n 


-MAX 

^<n 


are  best  bounds.  We  have  already  observed  that  QAD*  cannot  be  tricked  In  this  way,  so  these  may  not  be  best 
bounds  for  QAD*. 

The  worst  case  bahavlor  of  QAD  under  the  MIN  rule  Is  less  than  bout  half  »;  bad  as  under  the  AVE  and  MAX 
rules.  To  establish  this  we  will  first  need  two  lemmas. 


l**mma  4 - 1 .  Let  r^  denote  the  minimum  value  In  row  1  of  an  arbitrary  P  whose  rows  have  been  arranged  according  to 
the  MIN  rule.  Then 


*QAD  S  r1  +  2r2  +  ”•  +  nrn 


with  equality  only  If  the  r  occur  in  a  single  column  of  P. 

Proof.  Let  denote  the  contribution  of  Job  1  to  so  that 

^AD  “  C!  +  C2  +  •”  +  Cn’ 


If  we  esn  show  that  for  all  1,  Ct  £  ir^  then  (15)  will  follow  directly.  Thi  value  of  a  particular  C  will  be 
HJPiJ  f°r  ,0<”e  ^PP0**  th,t  rt  occur*  in  column  k,  so  ^  -  plfc.  Clear'.y,  hfc  £  1  and  rt  £  p^.  By  the  opera¬ 
tion  of  QAD  we  know  that  s  ,lfcPtk'  T*,en  **  ">*y  write 


C1  ■  hjpu  s  Vik  ■  Vi  s  lri 

and  (15)  Is  proved.  The  necessary  condition  for  equality  in  (15)  Is  obvious. 


Lemma  4-2.  If  r.  ir,  i  ...  lr  >0  then 
1  2  n 


r1  *  2r2  *  ...  *  nrn  ^  ^ 
rl  +  r2  +  ...  +  rn  2 


K 
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with  a quality  if  and  only  if  r  -  r  -  ...  -  r  . 

*  2  n 


£ro<?(-  Inaquallty  (16)  la  aqulvalant  to 


2+1  rl  +  2r2  +  ••»  +  nr  , 

2  '  r.  +  r,  +  ...■*  r  *  0 

1  2  n 


which,  alnca  all  ^  art  poaltlva,  la  ltaalf  aqulvalant  to 


(n+1)(r,  +  r2  ♦  ...  +  rn)  -  2(r,  ♦  2r2  +  ...  ♦  nrn>  *  0. 

Collact  tanas  In  (17)  to  gat 

(n-Dr,  +  (n-3)r2  +  (n-5)r3  +  ...  +  (5-n)rn  2  +  (3-n)ro1  +  (l-n)rn  s  0. 
Nov  match  r,  with  rB>  r2  with  r^,,  and  ao  on,  to  gat 


for  odd  n;  and 


(■•’)(',-'„)  +  <n-3)(VV,)  ♦  ...  ♦  (n-n)r ■-,*  0 

III 

(■•1)(rrrn)  +  +  ...  +  1(rn-rB  )  *  0 


2  T ' 


for  ....  to  to  con.tr.lnt  to.  r,  .  ...  ....  t.  „„  .„ 

(17)  to  th.r.for.  (IS),  .tolls,  Soil.  „,)  if  .„  0>  „ 

r,  -  r2  -  ...  .  rnj  otharvlaa,  .cm  tan.  will  ba  positive  and  (16)  will  ba  atrlct.  ■ 


Thao ram  4-2. 


_MIN 

-^2  <  211 

2  ' 


Proof.  Fir.t  notlca  that  If  tha  rf  ara  row  minima  of  P, 


F0PT  1  r1  +  r2  +  •••  +  rn 


.Its  .toll.,  ..1,  If  to  If.  1.  «ff .r.n.  column.  „  P.  (»„  u,to  „  n 

4-1  allow  ua  to  writ. 

*0*0  -  rl  *  2r2  +  *  nrn 

F0PT  rl  +  r2  +  -•  +  •  (21) 

Th.  oac.aa.ry  condition.  for  .quality  In  (20)  and  In  (15)  of  La-a  4-1  ara  contradictory.  Iha  ln.qu.llty  In  (21) 
1.  tharafo  a  atrlct.  U.lng  thl.  fact  tog.th.r  with  U—  4-2,  v.  g.t  (19)  dlractly.  ■ 

lh.  bound  in  (19)  ha.  not  ba.n  ahovn  to  b.  a  b.at  bound  for  alth.r  Q4D  or  QAD*.  In  fact,  tha  author  ha. 
baan  unabl.  to  conatruct  an  array  P  which  vlol.ta.  th.  aurpriain,  bound  in  tha  following  conjactira. 


Ihla  lamma  was  provad  by  Profassor  G.  W.  Stewart. 
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i  on  |ec  ture 

MIN 

K 

OPT 

We  now  turn  briefly  to  the  time  complexity  of  QAD  end  QAD*  with  row-sorting  rules,  tech  of  the  three  sort¬ 
ing  rules  examined  takes  time  proportionel  to  m  for  eech  of  the  n  rows  of  P,  or  O(ren),  to  celculate  the  sort 
keys.  Sorting  the  rows  tekes  en  edditionel  0(n  log  n)  ;  the  guts  of  QAD  require  O(rxn);  and  the  flnel  SPT  sort  of 
QAD*  takes  at  most  0(n  log  n) .  Both  elgorithms  with  sorting  rules,  therefore,  have  time  complexity 
0(max(mn,  n  log  n) ) .  Not  surprisingly,  this  is  the  seme  order  of  complexity  displeyed  by  the  mechlne  fsctor  al¬ 
gorithm.  For  purposes  of  comparison,  constructing  an  SPT  schedule  for  a  single  machine  -  perhaps  the  simplest 
scheduling  problem  of  ell  -  Itself  requires  time  of  0(n  log  n). 

5.  AVERAGE  PERFORMANCE 

Most  algorithms  beheve  much  better  most  of  the  time  then  they  do  In  the  worst  case,  end  the  QAD  elgorlthm 
Is  no  exception.  In  this  section  we  will  use  Monte  Carlo  methods  to  examine  the  everege  performence  of  QAD  end 
QAD*  under  verlous  row-sorting  rules  for  two  models  of  perellel  mechlne  systems.  Model  1  Is  extremely  simple 
and  correspondingly  remote  from  reellty;  Model  11,  e  computer  system  model,  is  considerably  more  down-to-eerth. 

Hie  meesure  of  performence  used,  es  In  Section  4,  Is  Fq^E/f0PT  (and  Fq^!/f0PT)  •  ln  e  typicel  experiment, 
a  large  number  of  errays  P  were  rendomly  genereted  e  id  the  values  of  this  performance  measure  calculated  for 
eech  of  four  sorting  rules:  MIN,  AVE,  MAX,  end  e  new  rule,  RAND,  which  errenges  the  rows  of  P  ln  ranoom  u-der. 
The  algorithm  of  Bruno,  Coffmen,  end  Sethi  [2,3]  wes  used  to  caJculete  FopT.  The  rules  MIN,  AVE,  etc.  used  with 
QAD*  will  be  denoted  by  MIN*,  AVE*,  utc. 

5 . 1  Model  1 

In  ,'hls  model,  the  p^  are  Integers  lnde-  .udently  drevn  from  e  uniform  distribution  over  a  specified  renge. 
Fig.  7  shows  the  results  of  en  experiment  in  which  were  genereted  200  errays  with  n  -  8  end  n  -  3,  where  the 
range  of  the  uniform  distribution  w< a  from  I  to  100.  Each  of  the  four  rules  was  epplled  with  QAD  to  eech  of  the 
200  problems,  end  the  four  curves  of  Fig.  7  ere  the  semple  cumuletlve  distribution  functions  for  the  velues  of 
FQADE/FOPT'  The  P°lnt  °*9°)  °n  the  AVE  curve,  for  example,  meens  that  90  percent  (180)  of  the  200  velues 

of  FQAD^FOPT  were  Fe,s  then  or  «9u*l  to  '-2.  We  may  Immediately  observe  thet  MIN  Is  substentlelly  better  than 
the  other  rules  for  this  experiment,  end  thet  the  performence  of  QAD  under  the  MIN  rule  is  quite  good:  39.3 
percent  of  the  semple  problems  were  scheduled  optlmelly,  and  95  percent  of  the  solutions  were  et  most  15  percent 
worse  than  the  optimal  schedule. 

The  performance  of  MIN  Is  little  Improved  by  the  flnel  SPT  sort  of  the  QAD*  algorithm,  though  the  other 
rules  are  Improved  to  'trying  degrees.  Fig.  8  shows  the  results  of  applying  QAD*  to  the  same  200  problems.  The 
curve  for  MIN*  la  omitted  because  Its  closeness  to  AVE*  would  heve  obscured  both  curves,  end  beceuse  It  Is  neerly 
Identical  to  MIN  in  Fig.  7. 

Table  1  shows  the  extent  to  which  each  rule  outperformed  each  of  the  others  in  this  experiment.  For 
Instence,  AVE*  yielded  e  better  schedule  then  MAX  ln  85  percent  (170)  of  the  200  sample  problems.  Table  2  gives 


MM 


-18- 


four  u..ful  ch.r.ct.ri.tic.  of  e.ch  ..ml,  di.tribution:  how  often  optimal  .ch.dul.  we.  found;  the  me.n  of 

the  dietribution;  the  v.lu.  et  end  below  which  95  percent  of  the  .en.pl.,  li.;  .nd  the  m.ximum  v.lu.  ob.erv.d  in 
the  experiment. 


TABLE  1 

Rule  A  better  then  Rule  B  (percentage) 


A\ 

B  MIN 

MIN* 

AVE 

AVE* 

MAX 

MAX* 

RAND 

RAND* 

MIN 

0 

0 

70 

30.5 

90.5 

51.5 

95.5 

55.5 

MIN* 

1 

0 

70 

30.5 

90.5 

51.5 

95.5 

55.5 

AVE 

19.5 

19.5 

0 

0 

82 

39 

86 

41.5 

AVE* 

27.5 

27 

73 

0 

93.5 

49 

94.5 

56.5 

MAX 

7 

7 

9 

5 

0 

0 

64.5 

21.5 

MAX* 

23.5 

23.5 

54 

24.5 

93.5 

0 

87.5 

43.5 

RAND 

3.5 

3.5 

12.5 

4.5 

33 

10.5 

0 

0 

RAND* 

21 

21 

52 

22.5 

76.5 

33 

96.5 

0 

TABLE  2 

Dietribution  Characteristic* 


PERCENTAGE 

OPTIMAL 

MEAN 

95 <  LEVEL 

OBSERVED 

MAXIMUM 

MIN 

39.5 

1.038315 

1.141304 

1.234177 

MIN* 

43 

1.038233 

1.141304 

1.234177 

AVE 

9 

1.087132 

1.228572 

1.516129 

AVE* 

38 

1.043284 

1.158228 

1.447581 

MAX 

1 

1.204095 

1.45946 

2.02963 

MAX* 

29 

1.069569 

1.220109 

1.357934 

RAND 

1.5 

1.280514 

1.614815 

1.917808 

RAND* 

26 

1.086  2  78 

1.259434 

1.449367 

The  evidence  thue  fer  preeented  euggeete  thet  for  Model  I, 
merglnelly  better  would  eeem  to  be  outweighed  by  the  eduitional 

i.tic.  of  the  result*  of  thle  experiment  turn  out  to  be  true  for 
thtt  will  be  described; 


MIN  1.  the  beet  rule.  The  feet  thet  MIN*  l, 
compu te tionel  effort  it  requires.  TVo  cherecter- 
•veri  experiment,  both  in  Model  I  and  Model  II, 


TO.  na..  U.W  1.  «...  of  ,0*0...  of  wrforu.c,  .„  HIM,  m,  Ktt,  UMD  ,„r 
2.  Hun  ~Cp.rf.nud  HIM  .im.ll,  MW*  ~.p.,f.n«d  MM»  .trtu.U,  th.  „„„ 

were  in  between. 

Severe  1  paremet.r.  of  th.  experiment  were  varied  to  test  th.  .en.itivity  of  th.  performance  of  QAD  with  MIN 
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”  Cl"*“  ■0d*‘-  *”  ""  ™.  ™  «.h  •  -  6  .«  -1.  ..h.r  „  ln 

fir.,  .„d  „ng.  of  dl.trIbutlon  „  ,.100„  jnd 

*0"’  *«  ch.r.tl.rl.tf., 

,h°"  »'  f‘!-  ’  ™“"*  -»  ,to  .pd 

(r.icrictions  ,.po..d  „y  tto  th.  tt>_ 

““  ‘"H—  '  -  ">■'-»  *»•.«..»  .»  ...... fo™  ,h... 

o.,.,  .d  „rf0„.„c,  m  ln  Ihrt<  >sptti_ 

,h*  -  -I.  -H  o,  ....  „  .  by  .  .  4. 

”“h  20  ”  ■"*  ”  Si  «“  »"k  . . .  .<  ....  6.  to  6.  *.  th. 

<"'K  *'*  «“  •-  »'»  o'  ««•  >  -  to  ..P.—  torto.„„.  ..... 

MS.  9.  „  .to,  ..  .  .... - -  -  -  ...  „>,  ,h.  ......  of  ,„«  dl„tlbulIo„.  „Mln.  nu 

stationary,  while  the  sample  variance  decreas^a  ivfl 

decreases.  The  three  means  are:  8  by  3,  1.038315;  16  by  4,  1.030334-  and 

>9  to  >.  1.0371205.  fto.h.r  by  ^  pf  /la 

“  6  ■•«».  •••»».  «  1.0303,9.  0PT 

difficult  to  conceive  of  a  reel  1. tic  situation  in  which  the  ere  Independent  and  uniformly  dis- 
A  job  which  1.  very  fast  on  some  machine  is  likely  to  be  fast  on  some  other  machine.  as  well;  a  , ma¬ 
chine  which  run.  one  Job  faster  than  the  other  machine,  is  likely  to  run  other  Job.  quickly  too.  We  turn  now 
to  a  more  realistic  model. 


5.2  Model  II 

Model  II  is  a  model  of  a  multiprocessor  computer  system  not  unlike  Crnegie-Mel Ion  University's  C..p  [„ 
^  m‘Chlne‘  are  dlfferent  m0del'  °f  ^  Dl8ltal  E’U1—  Corporation  PDP-11  computer  [6,.  and  the  Job.  are 

tors.  A  job's  m  proces.in*  times  are  calculate  by  dr.wing  a  sin,!,  time  from  an  exponential  distribution  and 

mult lplylng  by  the  machine  factors  associated  with  it.  inh  ri«  , _ . 

with  it.  Job  Class.  (Times  are  rounded  to  the  nearest  integer  ) 

lb.  Ml.  Met.  ««d  tor.  ....  c.lcul.t.d  „„  m.„  ,lg„„  (1) 

Hr.,  «Mt.-„  ,o  to  d..„,b«d  ,o..ld„.d  20.  prdbl...  ,uh  .to 

'  C“"  2’  J”b''  »'  “b‘'h  >1  -  «...  3,  .  ......  Jbb  .hloh  „old  „„ 

““  ’,°d'1  4“  b,"U"  <,bt  *”'P1,>  "  b*’"‘"d  *  ««»*•«  "-vie.  e«„.d  o.,y  ,0  ,to. 

Below  are  the  machine  factors  for  these  Job  classes: 


PDP-1 1/20 

PDP-1 l/4o 

PDP-1 l/4 5 

Class  1 

1 

.556 

.556 

Class  2 

1 

.291 

.134 

Class  3 

OO 

.556 

OO 

Note  that  floating-point  jobs  exploit 
greater  extent  than  "average"  Jobs. 


the  specialized  hardware  of  the  11/40  and  especially  the  ll/«  to  a  much 
The  expc  nflal  distribution  used  in  this  experiment  had  mean  1000 
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F"  “  ,h‘  -  — «•  «  p— .  «... ,  Job 

were  con.idered  to  be  it,  processing  time  on  the  1l/40. 

Th.  reeult,  of  thl,  experiment  ere  shown  In  Figure.  ,0  (QAD)  end  „  (QAD*)  end  In  Tables  3  and  4.  Comparing 

*’  “lth  Fl8'  7  W‘  °b,erV*  '  ^  Wlth  8Urp— *»*«  Performance  Is  better  In  Model  „  than 

-  Model  r  for  ell  rule,  except  RAND.  The  MIN  rule,  In  particular,  found  optimal  schemes  l„  over  60  percent 


A\ 

B  MIN 

MIN 

0 

MIN* 

0 

AVE 

13.5 

AVE* 

13.5 

MAX 

10 

MAX* 

11 

RAND 

0 

RAND* 

5 

TABLE  3 

Rule  A  better  then  Rule  B  (percentage) 


MIN* 

AVE 

AVE* 

MAX 

MAX* 

0 

48 

47.5 

71.5 

66.5 

0 

4d 

47.5 

71.5 

66.5 

13.5 

0 

0 

45.5 

37 

13.5 

18.5 

0 

56 

37.5 

10 

0 

0 

0 

0 

11 

11 

0 

47.5 

0 

0 

3 

2.5 

5 

4 

5 

15.5 

13 

32.5 

22.5 

RAND 

RAND* 

98.5 

89 

98.5 

89 

95.5 

79.5 

96.5 

81.5 

93.5 

63 

94.5 

70.5 

0 

0 

97.5 

0 

TABLE  4 

Distribution  Characteristics 


PERCENTAGE 

OPTIMAL 

MEAN 

95 i  LEVEL 

OBSERVED 

MA  VTMTTU 

MIN 

61.5 

1.01767 

1.081088 

riAALMUM 

1.303683 

MIN* 

61.5 

1.01767 

1.081088 

1.303683 

AVE 

39 

1.044154 

1.143041 

1.303683 

AVE* 

40.5 

1.039521 

1.139162 

1.303683 

MAX 

18 

1.088859 

1.257068 

1.303683 

MAX* 

24.5 

1.06943 

1.1884  95 

1.303683 

RAND 

1.5 

1.476827 

1. 952063 

2.249682 

RAND* 

4.5 

1.153185 

1.442467 

1 . 1?R7 

IKANirn 


i - 1 - 1 - 1 - r 

.1  1.2  1.3  1.4  1.5  1.6 


t 

Figure  1 1 


The  sensitivity  of  these  results  to  chenges  In  the  model' u  perimeters  was  examined  by  running  several  addi¬ 
tional  experiments.  As  in  Model  I,  changing  the  range  of  values  of  the  p^  had  little  effect.  An  experiment  of 
100  sample  arrays  was  run  with  all  parameters  as  before,  except  that  the  mean  of  the  exponential  distribution  was 
changed  from  1000  to  60000  (milliseconds).  The  performance  of  QAD  under  the  MIN  rule  was  not  significantly  dif¬ 
ferent  from  that  of  the  first  experiment.  Again,  strong  conclusions  are  ruled  out  by  the  limited  extent  of  the 
testing,  but  performance  In  this  model  would  tentatively  appear  not  to  depend  on  the  exponential  distribution. 

In  further  agreement  with  Model  I,  performance  was  found  to  be  sensitive  to  changes  In  n.  An  experiment  was 
performed  In  which  the  number  of  Jobs  In  each  class  was  doubled,  bringing  the  total  to  16.  A  fourth  machine  was 
added  to  the  model,  a  PDP-11  Model  05,  whose  machine  factors  for  the  three  Job  classes  are:  Class  1,  1.25;  Class 
2,  1.25;  and  Class  3,  “.  One  hundred  sample  arrays  were  -ested.  The  results  appear  In  comparison  with  the  8  by 
3  results  In  Fig.  12  on  a  greatly  expanded  horlaontal  scale.  The  means  of  the  two  experiments  are:  8  by  3, 
1.01767;  16  by  4,  1.019837.  Approximately  the  same  behavior  found  In  Model  I  Is  demonstrated  here:  as  n  In¬ 
creases  (with  m  ■  log  n),  the  sample  mean  appears  to  remain  relatively  stationary,  while  the  sample  variance  de¬ 
creases. 


Observed  Prob  (F 


Figure  12 
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6.  CONCLUSION 

For  the  scheduling  problem  studied  In  this  paper,  the  QAD  algorithm  seems  to  be  a  reasonable  alternative  to 

the  optimal  algorithm  of  Bruno,  Coffman,  and  Sethi  [2,3].  The  QAD  algorithm  (with  the  MIN  rule)  takes  computa- 

2  3 

tional  time  of  0(max(mn,  n  log  n) ) ,  while  the  optimal  a  lgoritlm  requires  0(max(mn  ,  n  )).  Furthermore,  QAD  Is  an 
extremely  simple  algorithm,  and  easy  to  work  by  hand  (as  might  be  required  In  an  industrial  shop,  for  Instance). 
By  contrast,  the  optimal  algorithm  Is  quite  difficult.  While  QAD  does  not  always  find  optimal  schedules,  It 
frequently  does  (and  It  always  can) ;  and  Its  performance,  bounded  In  the  worst  case,  appears  from  limited  experi¬ 
mental  results  to  be  very  good  most  of  the  time. 

There  Is,  of  course,  much  room  for  further  work.  The  most  Interesting  unanswered  question  Is  whether  here 

exists  a  simple  row-sorting  rule  which  will  guarantee  QAD' s  production  of  optimal  schedules.  It  might  be  the 
case,  however,  that  sorting  according  to  any  function  of  a  row  Is  by  Itself  insufficient;  more  Information  about 
the  processing  time  array  might  be  required  to  discover  the  optimal  row- permutation  promised  by  Theorem  3-1.  If 
im  optimal  rule  cannot  be  found,  perhapa  row-sorting  rules  more  fruitful  than  the  simple  ones  considered  here 
cxild  be  discovered.  A  proof  of  the  conjecture  In  Section  4  would  be  extremely  Interesting.  Finally,  the  com¬ 
putational  complexity  of  this  problem  is  unknown,  and  the  work  reported  here  only  begins  to  suggest  that  the  com- 
3 

plexity  Is  less  than  0(n  ). 
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